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^ ' A simple model of a two-mode non-resonant parametric amplifier is studied 

c3 I with special regard to non-classical features such as revivals and squeezing. 

The methods used apply for an arbitrary pump parameter. Detailed analytical 
and explicit expressions are given when the coupling of the two modes has an 
harmonic time-dependence. Despite its simplicity the model exhibits a very 
broad range of intricate physical effects. We show that quantum revivals are 
possible for a broad continuous range of physical parameters in the case of 
initial Fock states. For coherent states we find that such revivals are possible 
only for certain discrete rational number combinations of the ratio of frequency 
detuning and pump parameters. Correlation effects are shown to be very 
sensitive to the initial state of the system. 
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1 Introduction 



The quantum-mechanical analysis of coupled harmonic oscillators enters into many 
fields of physics. In the design of detectors for gravitational radiation the dynamics of 
an harmonic oscillator and a transducer is e.g. of importance (see e.g. Ref. 0). The 
coupling of free quantum fields to classical external fields can also be thought of as a 
system of harmonic oscillators with, in general, space- and time-dependent couplings (see 
e.g. Ref. [0]). In the field of quantum optics similar systems of coupled harmonic oscillator 
systems are of great interest as e.g. in the description of parametric down-conversion of 
photons (see e.g. Refs.g |, |, |, 0, |, |, 0, [TT], |T1). 

Our primary goal here is to study a simple interaction between a LC-circuit and one 
mode of the second quantized radiation field all without damping effects included. The 
corresponding free Hamiltonians are used to define reference states as e.g. the notion of 
a photon. In order to calculate photon-number distributions at any time normal ordering 
procedures become an essential ingredient in our study. In particular we have, in terms 
of a canonical rescaling of variables, considered the following classical Hamiltonian 

H = \PI + \ulQl + \P^ + \ulQl - g {.[^QaP, + .[^PaQl 
Z Z, Z, Z, \y UJi) y UJa 

where the canonical momentum Pa is proportional to the electric current / = dQa/dt of 
the LC-circuit with inductance L and capacitance C {uJa = 1/yLC). Pb and Qi, describe 
the harmonic oscillator of the singe-mode of the radiation field with frequency uoi,. The 
last interaction term is a simple form for the interaction of the radiation-field with the 
LC-circuit in terms of a real coupling constant g. The particular form of the interaction 
term in Eq. ([1|) is motivated by the fact that a minor extension of this model can also 
be used to describe non-resonant parametric down-conversion, as will be discussed in 
more detail Section ^ Eq. (|lD is also the Hamiltonian for a two-dimensional electrically 
charged harmonic oscillator, in general non-isotropic, with an external magnetic field 
B = g{^Jujb/uJa — \/^a/^b) perpendicular to the ab-plane and with angular frequencies 
cj^ — > c<j^ — g^uJbjijJa and uJl ^ ujI — g^uJajijJh- If the original angular frequencies uJa and lo^ 
are equal the harmonic oscillator becomes isotropic and the magnetic field is zero. 

Despite the very simple structure of the Hamiltonian Eq. (Q) we will see that the 
corresponding quantum dynamics can lead to quite intriguing physics. In view of the 
impressive experimental developments in quantum optics the system under consideration 
may actually be realized in the laboratory. 



The paper is organized as follows. In Section 2 the dynamical system is defined in more 
detail. An exact treatment of the time-evolution operator is presented in Section 3. The 
results obtained turn out to have a broad range of applicability. Quantum revivals and 
their dependence of initial states are discussed in Section 4. In Section 5 non-classical 
features are studied in terms of correlation functions. The effect of detuning on the 
squeezing properties of a parametric amplifier is studied in Section 6. Signal-to-noise 
aspects are discussed i Section 7 and some final remarks are given in Section 8. 

2 The Dynamical System 



The coupled harmonic oscillators we consider are described by the two-mode Hamil- 
tonian 

H = Ha + Hb + Hint , (2) 

where 

Ha = UJaiUa + -) , 
Hb = UbiUb + -) , (3) 

in natural units {h = c = 1) and where n^ = a'^a and Ub = b'^b. The interaction between 
the a and b mode is chosen to be of the form 

H,nt = ^[9it)ab-g*it)a^b^] , (4) 

where g{t) is a time-dependent coupling. If g(t) = g we recover the system described 
by Eq. (|1|). The interaction Hamiltonian Eq. (^) describes a two-mode non-degenerate 
parametric down-conversion process, where the modes are known as "signal" and "idler" 
and where the pump parameter g{t) then describes an arbitrary classical pump field (for 



reviews see e.g. Refs.y, 0, [T|, |T|, |T§). 



The Hamiltonian (^ can be written in the form 

Hi^t = i [g{t)K_-g*it)K^] , (5) 

where 

K^^ab , K+^ a^b^ , K^ ^ 1^±^ , (6) 



span the Lie algebra of the non-compact group SU{1, 1), i.e. 

[K+,K^] = ~2Ko , [Ko,Ki] = ±K± . (7) 

The Casimir operator of the Lie-algebra Eq. (|^) is given by C = K^ — K^K_ — Kq = 
$($ + 1), where we can choose —2$ =| Ua — nf, \ +1. We therefore see that we have a 
realization of the unitary representation D+($) of SU{1, 1) [jl^. We observe that when 
we make the identification a = b, the Hamiltonian describes a degenerate parametric 
amplifier. This corresponds to a realization of the SU{1, 1) Lie-algebra Eq. (^ with 
K_ = a^/2, K^ = (a^)^/2 and Kq = {a'^a + aa^)/4 which also constitutes a realization of 
a discrete unitary representation of 5'^(1, 1) ||T^. 



3 The Time-Evolution Operator 

The Hamiltonian in the interaction picture and in terms of the operators K^ and K_ 
reads 

H^^^{t) = i [^(t)e-^(^"+'"'')*ir_ - ^*(t)e*('^"+'^'')*K+] , (8) 

with the corresponding Schrodinger equation 

^JUI{t) = t [ ~g{t)K- - r{t)K+ ] Uj{t) , (9) 

and where we have defined g(t) = g{t)e^''^'^'^~^'^''^^. An exact analytical expression for the 
time-evolution operator Uj(t) can be obtained as follows. Due to the presence of the 
Lie-algebra of SU{1, 1) we employ the Wei-Norman technique |]TB|, |T3 and express the 
time-evolution operator in the following form 

since this is a convenient choice when computing expectation values and probability ampli- 
tudes. If we substitute Eq. ( p!OD into the Schrodinger equation Eq. (H) we find that A^{t), 
A_{t) and Ao(t) must satisfy the following system of non-linear differential equations: 

^^A4t) = mAi{t)-r{t), 

j^Aoit) = mA^t), (11) 
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with the initial conditions A+(0) = A-(0) = Aq{Q) = 0. These initial conditions ensure 
that Ui{0) = 1. Eqs. (|TT1) are actually valid for the Hamiltonian Eq.(^ in general and 
can e.g. be studied by making use of various numerical methods. In the present paper we 
will, however, consider a simple situation in which case analytical methods are at hand, 
i.e. we consider a coupling g{t) with harmonic time-dependence 

9{t) = ge'^' , (12) 

where g is a real constant. For parametric down-conversion processes the frequency u is 
then referred to as the pump frequency. In the course of preparation of the present paper 
we have observed that in the special case of the time-dependence of Eq. (|12D algebraic 
methods can be used to solve for Eq. (||) (see e.g. Ref . [^ ) . To the extent one can compare 
our methods and results with those of Ref.[^] we find agreement. We, however, once more 
stress that our methods in principle can be used for an arbitrary time-dependence of g{t). 
It is convenient to define 

Q = UJ — Ua — UJb . (13) 

When the frequencies Ua and ui, sum up to the classical pump frequency u, i.e. i7 = 0, 
we have parametric resonance. In this case the Schrodinger equation is easy to solve (see 



e.g. Ref. 11). 

Let us now consider the more general situation i7 7^ 0, in which case the interaction 
Hamiltonian constitutes a simple model of a detuned "broad-band" pump in parametric 
amplification. Within first-order perturbation theory, in which case the time-ordering 
procedure of the time-evolution operator is not important, the broad-band pump down- 



conversion has been considered in great detail in the literature ||2l], ^ |^, ^, ^ |2^, ^, ^ 



2U]. Here we are interested in an exact treatment of the time-evolution operator, i.e. we 
have to solve Eqs. (|TT]) with the harmonic time-dependence Eq. (|T^. An exact treatment 
can, e.g., be used to investigate unitarity effects (i.e. including not only the vacuum 
and the two-photon states but all possible intermediate states) of photon correlation 
experiments. 

A important point here is that we have to solve Eqs. (0) separately for three different 
cases: fc^ < 1, fc^ > 1 and fc^ = 1, where we have defined the dimensionless parameter 

k^^ . (14) 



In the case fc^ < 1, the solution has the form 



A At) 



— e 



Ant 



Vl — fc^ tanh ( gt^/l — k^ — ^7 ) + ik 



A_(t) = Vl^^ tanh (gtVl^^ ~ i'j) + ik , 



(15) 



Ao{t) 



-In 



cosh {gWl — A;2 — 27) 



cos 7 



iQ 

—t 

2 



where 



tan 7 



k 



(16) 



The solution in the case fc^ > 1 is of the same form as in the case above with the 
substitutions 7 — > — 7(= 5), tanh — »• cot and cos — >• ?sinh. The explicit solution of the 
Eqs. ( p^D is then 



A At) 



— e 



'iUt 



VA;2 - 1 cot (gWk'^ - 1 + 'i5 ) + zfc 



A_(t) = Vfc2 _ 1 cot (gWk^ -l + i6) +ik 
sin {gty/W^^ + i6) 



:i7) 



^(t) 



In 



sin(i5) 



in 

—t 
2 



where now 



coth5 



VP 



;i8) 



If fc^ = 1 it is particularly easy to show that 



A At) = - 


_ ^-int 


1 
gt + i 


A At) = 


1 

:+i 

_gt + i 


) 


Aoit) = - 


- In [1 - 


ig\ 


^ 2 



(19) 



is the solution of the Eqs. (|TT]) with the correct initial conditions. We notice that A At) = 
— e~*^*A_(t) for all three cases. 

Due to the constraint of unitarity, i.e. UAt) = UAit), the functions A^(t), A^{t) 
and Ao{t) are not independent. By making use of operator reordering techniques (for a 



pedagogical account see e.g. Ref.[Q), the unitarity of f//(t) leads to the general conditions 



(20) 



r^j^\Lj — 


e2AoW_^_(t)A+(t) 


A*_{t) = 


A^{t) 


e^Aoit) ^ A_{t)A+{t) 


Ao{t)-A*{t) ^ 


l-A_(t)A+(t)e-2^°W 


Ao{t)-A^{t) ^ 


A- it) A^t) 



and hence 

Ao{t)-A*g{t) ^ _ _ 

A;(t) A*4t) ■ 

Using (1^) the following useful relationship follows 



(21) 



-Ao(t)-A^,(t) [l_|^_(t)|2] =1 . (22) 



It is not obvious that the explicit solutions Eq. (|T5|) - Eq. (|1^) for A^(t), A^(t) and Ao(t) 
satisfy the Eqs. (EH). However, by using the fact that A^(t) = — e~*'^*y4_(t) and various 



trigonometric relations, one can explicitly show that Eqs. (pO]) are fulfilled for all three 
cases k"^ < 1, k"^ > 1 and k"^ = 1. 

With the solutions above it is now straightforward to calculate various expectation 
values and transition amplitudes. In this context it is useful to obtain the annihilation 
operator anit) = a(t) in the Heisenberg picture (a(t = 0) = a), i.e. 

a{t) = U{tyaU{t) , (23) 

where 

U{t) = Uo{t)Uj{t) , (24) 

and 

Uoit) = exp[-i{ujaa^a + ujbb^b + Eo/2)t ] . (25) 

Here Eq = Ua + ujb and f//(t) is given by Eq. (|10|). Straightforward calculations now gives 

a{t) = e-^5-*^''* [a - b^ A*_{t)] , (26) 

from which b{t) follows by the substitution a ^^ b. The creation operator a^ follows from 
Eq. (p6D by taking the Hermitian conjugate. When k'^ = this solution for the operator 
a(t) reduces to the wellknown result (see e.g. Ref . p^ ) 

a{t) = acosh{gt)-b^smh{gt) . (27) 



The commutator [a{t),a^{t)] is in general given by 

[a(t),at(t)] = e-^oW-^SW [i _ |A_(t)p] . (28) 

Due to Eq. ( p2D the quantum condition [a(t), a"''(t)] = 1 is satisfied for all times, as it 
should be. 

The Hamiltonian of the system is such that a^(t)a(t) — 6"''(t)6(t) is a conserved quantity. 
We obtain 

Uait) = a\t)a{t) = 
a^a - A_{t)K_ + A+(t)e-2^«W {K+ - 2A_{t)Ko + Al{t)K_) . (29) 

Since nt,{t) = b'^{t)b{t) can be obtained from Eq. ( p9|) by means of the substitution a *-^ h 
it follows immediately that riait) — rihit) = na{0) — "^,(,(0). The time-dependence of other 
combinations of annihilation and creation operators can be obtained in s similar manner. 



4 Quantum Revivals 



In this section we focus our attention on quantum revivals, i.e. when the initial state 
is reproduced by the time-evolution (for a reviews see e.g. Refs.|l3l[|). The time-evolution 
operator Eq. (|lO|) is written in a form that makes the evaluation of transition amplitudes 
(/I U{t) \i) straightforward. It is clear from the Eqs.(|T^ that exact quantum revivals will 
occur for Fock states only if fc^ > 1 and at revival times t^ev such that 



gtrevVk"^ — 1 = nvr , (30) 

where n is a positive integer. This is so since the corresponding transition probability is 
independent of Q. The phase in the transition amplitudes due to the free time-evolution 
will not effect the transition probabilities, at least not for initial and final Fock-states. For 
coherent states we will, however, find it convenient to evaluate the transition probabilities 
in a basis in which case the free time-evolution is absent. For coherent states it then turns 
out that Eq. (|30D is extended in such a way that only certain rational numbers of k"^ will 



lead to revivals. Here we therefore see a remarkable difference between initial and final 
Fock states and coherent states respectively with regard to quantum revivals. We now 
treat these special cases in more detail. 



4.1 Fock States 

The two-mode initial state has now the form \i) = \r)a ® \s)f, = \r,s). We also use 
the notation (/| = ;,(r7i| ^ a{n\ = {m,n\. The transition amplitudes to be derived are 
Cmn(t;r,s) = {m,n\Ui{t)\r,s). The probabilities to be calculated are then Pmnit;r,s) = 
\cmn(t;r,s)\'^. We will suppress the initial state indices r and s if it is clear from the 
context what initial state we are considering. The corresponding amplitudes Cmnit) can 
be found in closed form by making use of the results from Section ^, i.e. 

^ {r-k)\{s-k)\ k\ {n + k-r)\ ' ^ ' 

It follows from the actual expressions for Ao(t), A+(t) and A^{t) that |cmn(i^)P does not 
depend on Vt and therefore one obtain the revival times Eq. (|30D . 

We now present more explicit results in a few simple cases. If the initial state is the 
vacuum state, i.e. r = s = 0, it follows that 

P^n{t) = Smn e^°W+^S « (| A+(t) H'^ . (32) 

It can now be shown explicitly that this probability distribution is properly normalized, 
i-^- Ylmn=oPmn(t) = 1 for all times. Pmnif) can be explicitly evaluated by making use of 
trigonometric identities. One can then derive the useful results 

(33) 
(34) 

(35) 



and 



a 


= A. 


e-Mty 
.it) ' = 
-- -A (1 


-Alit) ^ 

\A.{t) ' 
'\A. (f\f 


sinh^(c/tVl -k'^) + l 


-P 


) 


y{t) 


2 sinh.'^ {gt\/ 1 — 


P) 




sinh^(^tVl - A;2) + 1 - 


A;2 


n 


,-2Ao(i) _ sinh"(^tVl - 


-F) 





valid for all k"^ . The reason why we use the notation no(t) is the fact that Eq. (|35D is noth- 
ing else but the expectation value of the number operator na{t), as given by Eq. (^) (or 
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nb{t)), in the initial vacuum state. We notice that the last three equations are connected 
in the following way: x{t)y{t) = no(t) and x(t) = 1 + no{t). 

The reduced density matrix for the a-modes (or the 6-modes) is again given in terms 
of Eq. (|32D. As has been noticed before (see e.g. Refs.p, |^, ||, ^, 0) the reduced density 
matrices describe a system at finite temperature T = l/f] such that exp(— /^cj) = y(t), 
where u = Ua (or u = Uf, for the 6-modes), i.e. if fc^ < 1 the temperature approaches 
infinity in the limit of large gt. This interpretation has, of course, a limited range of 
validity. 

With \i) = |1, 1), the corresponding probabilities are given by 



P„„W = «™„|j^y"-W-^."-'W + y"«WJ;;(^ , (36) 

and again it follows that Pmn{t) is properly normalized. Furthermore we have {i\na{t)\i) = 
1 + 2no(t). The probabilities in Eq. ( p6D for n = m = 1 and n = m = 3 are shown in 
Fig. |l| for k"^ = 1.5 in which case quantum revivals are predicted. In Fig. ^ the same 
probabilities are shown, but now for the case k'^ = 0.5. These figures reflect the typical 
physical situation: For k"^ < 1 the probability has initially a peak structure and then it 
decays exponentially with time. Actually all probabilities Pmn{t) will then tend to zero for 
sufficiently large times in such a way that their sum still is one, i.e. J2m n=o P"in(t) = 1- On 
the other hand, if A;^ > 1 the probability exhibit a totally different oscillatory behavior for 
all times with quantum revivals as predicted by Eq. (PUJ). The reduced density matrices 
can be computed as above. For large gt and A;^ < 1 one finds that the reduced density 
matrices correspond to the same thermal distribution as for the initial vacuum state as 
discussed above. 

Let us now assume that there is initially no 6-quanta, that is \i) = \ip)a ® |0)fe, where 
the initial a-mode state is a general pure state described by 

oo 

\i^)a = Y,VVse"^^\s)a . (37) 

Here Vg is the probability that the a-mode is in the Fock state \s)a and (ps is the corre- 
sponding phase. As above, the final states Uo\r,s) and \r,s) give the same result for the 
probability. The probability Pmn(^) = \{Tn,n\ Ui{t) \iIj,0)\'^ is given by 

^""^^ "-"^(n-m)!m!x("-™+i)(t) ' ^ ^ 

independent of the phase 0„. Again one can analytically verify that J2m n=o P"iri(t) = 1 for 
all times when ^^q Vn = ^ for all three cases k"^ > 1, k"^ < 1 and k"^ = 1. The expectation 



Pmn\J' ) 




Figure 1: The probabilities pii{t; 1, 1) (solid curve) and P33(t; 1, 1) (dotted curve) 
when k'^ = 1.5 and with an initial vacuum state. The corresponding revival-times 
are integer multiples of 7r\/2 ~ 4.44. 



values of the number operators for the a and 6-modes are {i\na(t) \i) = {ua) +no{t) [{ua) + 1] 
and {i\nf,(t)\i) = nQ{t)[{na) + 1] respectively. Here we have that (n^) = J2'^=o^'^n- If in 
particular the a- mode is Poisson distributed the probability Pi2(t) is shown in Fig. ^ 
for /c^ = 0.5 and k"^ = 1.5. When k'^ > 1 the probability oscillates between and 
the maximum value oi pi2{t), i.e. 8|apexp(— |ap)/27 which is close to 0.104 for the 
parameters corresponding to Fig. |^. Furthermore, the case k"^ = 1 is very much like case 
k'<l. 

The reduced density matrices for the a- and 6-modes are again diagonal with the 
matrix elements 



Pm{b;t) = ^Pmn(t) = X] '^' 



n=0 



1=0 



m! /! x'+^t) 



(39) 



for the 6-modes and 



Pn{a;t) = ^Pmn{t) 



m=0 



X 



n+1 



W; 



m=0 




<it) 



(40) 
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PmnyJ' ) 




gt 



Figure 2: The probabilities pu{t; 1, 1) (solid curve) and P33{t; 1, 1) (dotted curve) 
when k"^ = 0.5 and with an initial vacuum state. 

for the a-modes. For large gt and ii k^ < 1 one finds again that both modes are thermally 
distributed with the same distribution as in the case of an initial vacuum state. 

4.2 Coherent State 



We now consider initial and final coherent states. The coherent state transition 
amplitudes are somewhat general in the sense that all Fock-state amplitudes can be 
obtained from them in a standard manner (see e.g. Ref. ||32||). If the final state is 
a coherent state, then the presence of Uo(t) is important since Uo(t) \(y)a ® \P)b = 
g-iBot/2 |«e~*'^"*)a® |/?e~*'^''*)b. To exhibit quantum revivals we therefore find it convenient 
to use Uo{t) \a)a ® \l3)b as a final state. Let us, in particular, consider the case when the 
system is initially in the coherent state \i) = \a)a ® \(3)b = |«,/5) and final state (/| = 
i,{z\ a{uj\ f/o(t) = {z,w\ t/o(t). The transition probability p^^, (t) = \{z,w\Uj{t)\a, /3)\'^ is 
then given by 
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gt 



Figure 3: The probability Pmn{t) in Eq. (pSD when n = 2 and m = 1 for the two 
cases k'^ = 1.5 (solid curve) and A;^ = 0.5 (dotted curve) and when the a-photons 
initially are Poisson distributed Vn = |ap"e~'"' /n! with amplitude a = 0.85. 



Pzw{t) 



^ . A+(t)w*;2*|2, A_(t)a/?|2|X«e"^''^*^2|X/?e^°^*^2 

/n O O O O 



(41) 



Eq. (^TD can be used as a generating function for the Fock state transition probabilities. 
The expectation value of the number operator for the a-mode is 

{i\na{t)\€) = {ria) + no{t)[{n,) + (n^) + 1] - 2Re[a(3A_{t)]{l + no{t)) , (42) 

and for the 6-modes {i\nb{t)\i) = {i\na(t)\i) + {ub) — {ua). Here we have {ua) = |ap 
and {rib) = |/5p- In contrast to the Fock case, this transition probabihty depends on the 
detuning frequency Q. To obtain exact quantum revivals, the parameter k"^ can therefore 
not take an arbitrary value. To investigate the revival period of Eq. ( ^T]) it is convenient 
to rewrite the revival probability in the following form 

g-(|ap + |/5|2)[2-e^oW-e^S(*)] 



PPait) 



X e 



aP[A4t) + A*^{t)]^2Mt) + Al{t) 



(43) 
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By inspection of the actual expressions for A+(t) and A^(t) in Eq. ( ]T7| ) we see that the 
second factor in Eq. ( ^3D is periodic when 

where n and p are integers satisfying n > p, i.e. k'^ has to be a rational number in order 
to obtain exact quantum revivals. From the first factor Eq. (^31) we also see that Ao(t) 
has to be periodic in order to observe revivals. Thus, we have further restrictions on the 
integers n and p: if n is even, then p must also be even. If n is odd, then p must also be 
odd. The revival time corresponding to Eq. (|^) is then 



gtrev = TT ^/ n? - p"^ . (45) 

The solid curve in Fig. ^ shows the probability p/j^ (^) when /c^ = 9/5, i.e. when k"^ satisfies 
Eq. (^4|) {a = P = 1). The revival times are then integer multiples of gtj-ev = 7rv^20 ~ 
14.05. 

When the parameter fc^ does not satisfy Eq. ( ^Ij ) the time- dependence of ppait) is 
much more complicated. As long as |a| 7^ and |/3| 7^ the probability in Eq. ( ^Bf ) will 
then never reach unity. Hence, there are no revivals. But as we can see from Fig. ^ the 
probability ppa{t) nevertheless exhibits some peaks at e.g. gt ~ 22 and gt ^ 41 (dotted 
curve). These times correspond to the condition 2 — exp[ylo(t)] — exp[74Q(t)] ^ 0, due to 
the first factor in Eq. (^). If this condition is not satisfied Ppait) will be exponentially 
suppressed when a and P are large, i.e. when the photon number is large. 



5 Non-Classical Behavior 

Our two-mode system exhibits quantum correlations. These correlations can be de- 
scribed by various correlations functions. We find it convenient to consider the following 
well-defined second-order correlations function: 

fit) = yWW)vWW) - {a^a b^ , (46) 

where e.g. 

(a^ab^b) = {i\ U\t) a^a b^b U{t) \i) . (47) 
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Pfioc{t) 




gt 

Figure 4: The solid curve shows the probability ppai^) when k'^ = 9/5 and 
a = j3 = 1. The dotted curve shows ppait) when /c^ = vr and a = (3 = 5. The revival 
times corresponding to the solid curve are integer multiples oigtrev = vrv20 ~ 14.05. 



Since the expectation-values in Eq. (^) are normal-ordered we can always write these 
moments in terms of a Glauber-Sudarshan P representation, e.g. 



{a^ah^h) 



S(x / S(5 \(x\ 



P(«,/3;t) 



(48) 



If P(a, /5; t) is real and positive, then the moments in Eq. (|46|) are defining inner products. 
Hence, the Cauchy-Schwarz inequality implies the /(t) > 0. States described by such a 
function admit a classical interpretation where the complex field amplitudes a and /? 
may be treated as stochastic random variables with the joint probability distribution 
P{a,l3]t). If, however, f(t) < then P{a,l3;t) is not real and positive everywhere. In 
such a situation the correlations cannot be described classically since P{a,P;t) cannot 
be interpreted as a classical probability distribution. Hence, the sign of f{t) reflects the 
classical or quantum nature of the correlations between the modes. 

Another convenient way to characterize a non-classical state for single modes is the 
Mandel quality factor Qa defined by ||15[ 



Qait) 



Mt)) 



(49) 
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where na{t) = a^{t)a{t) and (AC(t))2 = Var[0{t)] = {0\t)) - {0{t))^ for an observable 
O. The quahty factor is constructed in such a way that Qa = when the statistics 
describing the isolated a-quanta is Poissonian. If the Qa factor is negative the state is 
quantum mechanical, without any classical analog. In this case the statistics is sub- 
Poissonian. The smallest value —1 is obtained for a Fock state. On the other hand, when 
Qa is positive the statistics is super-Poissonian and the isolated a-mode can be described 
in classical terms. A corresponding formula as Eq. ( ^Of ) holds for the 6-quanta. 

5.1 Fock States 



Let us now consider the situation with initial Fock states \r)a ® \s)b = \r,s). The 
Mandel factor in this case is given by 

_ nojt) 2rs + ng(t) [2rs + r + s + 1] - r 
r + no(t) [r + s + IJ 

where no(t) is as in Eq. (|35|). If r = then (5a(0) = 0, independent of s, see e.g. Fig. ^ 
(dotted curve). On the other hand, if r 7^ then (5a(0) = —1, also independent of s, see 
e.g Fig. ^ (thin solid curve). The physics when t > depends strongly on the value of 
k"^. li k^ < 1 the Mandel factor increase exponentially with time, which leads to super- 
Poissonian statistics of the a-quanta. If, on the other hand, k"^ > 1 then the (5a(^)-function 
will have an oscillatory behavior. The maximum value of Qa{t), when r 7^ and s = 0, 
is Q™"^ = [1 + 1/r - (A;2 - 1)2]/[(P _ ^^2 ^ (^^2 _ ^^^^ ^ ^/^^j^ If r = and s ^ then 
Q^o.^ _ 1/(^/^2 _ ^y rjnj^g minimum value of Qa{t) is or —1, depending on the value of r 
as discussed above (c.f. Fig. ^). 

The /(t)-function when the system is initially in the Fock state \r, s) is given by 

fit) = x\t) ^r{r - 1) + Ar{s + l)y{t) + {s + l){s + 2)y\t) 
X ^/sis - 1) + 4s(r + l)y{t) + (r + l)(r + 2)y^{t) 
- x^{t) [rs + {r + l){s + l)y^{t) (51) 

+ {rs + r{r + 1) + s{s + 1) + (r + l)(s + 1)) y{t) ] , 

where x{t) is given by Eq. ( |5B| ) and y(t) = 1 — l/x{t) as in Eq. (0). Notice that f{t) is 
symmetric with respect to exchange of r and s as it should. 
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Fit) 



Qait) 




4 gt 6 



Figure 5: The Mandel factor Qa{t) when the system is in the Fock state \r)a'Si\s)b 
and when k^ = 1.5. The thin sohd curve correspond to r = 1, s = and the dotted 
curve correspond to r = 0, s = 1. The thick sohd curve shows the corresponding 
F(t)-function. 

Let us now, for reasons of convenience, consider the correlation function F(t) = 
f (t) / >y {na{t)) {ni,{t)) instead of f{t) itself. Fig. | shows F{t) as well as the correspond- 
ing quality factors for k"^ = 1.5 and r = 1, s = or r = 0, s = 1. The F(t)-function 
in this figure is negative for all times, except at the revival times where it is zero, and 
the two-mode down-converted light therefore exhibits quantum mechanical correlations 
in this case. Since the Fock state is a typical quantum mechanical state this result is 
as expected. Nevertheless, the quality factor oscillates between and 2 (dotted curve) 
indicating that the statistics of the isolated a-system is super-Poissonian. It is, however, 
possible to make F{t) positive in spite of the fact that the Fock state is a typical quantum 
mechanical state. This can be done by choosing r ^ s 7^ or vise versa (e.g. r = 50 and 
s = 10 as in Fig. |^ (thin solid curve)). If r = s then F{t) = — 1 for all possible values 
of k"^, c.f. the thick line in Fig. ^. If s = then F{t) < also for all fc^. Fig. ^ shows 
the typical behavior when k"^ < 1. As we can see, F{t) approaches —1 for sufficiently 
large times. Furthermore, the thin dotted curve in Fig. |^ is negative for all times, but the 
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F{t) 




gt 



Figure 6: The normalized function F{t) = f {t) / ^/njfynijf) when k'^ = 1.5. The 
thin sohd cure shows F{t) when the system is initiahy a Fock state \r)a ® |s)6 with 
r = 50, s = 10 and the dotted curve when r = 50, s = 0. The thick sohd curve 
corresponds to arbitrary r = s. 



thin solid curve has a positive peak at small times (0.1 ^ gt ^ 1.3). In this short period 
of time the correlation between the a- and 6-mode can be described in classical terms. 
For all other times, when F{t) is negative, the down-converted quanta exhibits quantum 
mechanical correlations. Furthermore, the corresponding Mandel factor starts at —1 or 
(depending on the value of r) and increase exponentially with time. 

As mentioned earlier, in the case fc^ < 1, the isolated a- or 6- mode describes a system 
with a temperature that approaches infinity in the large gt limit. It is therefore expected 
that the statistics become super- Poissonian in this limit. 



5.2 Coherent States 

Since the quantum mechanical state that corresponds to a classical description is the 
coherent state, we expect classical-like behavior when considering initial and final coherent 
states. Fig. |] shows the function F(t) for the initial coherent state \a)a <S) \l3)b = \c(,l3) 
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F{t) 



Qait) 




Figure 7: All curves correspond to k'^ = 0.5. The thin solid curve shows F{t) 
when the two-mode system is initially in the Fock state \r)a ^ \s)f, = \r,s) with 
r = 50, s = 10. The thick solid curve is the corresponding Qa{t)-ia,ctor. The thin 
dotted curve shows F{t) in the case r = 50, s = and the thick dotted curve is the 
corresponding Qa{t)-iactor. 

when k"^ = 1.5. For reasons of comparison we have chosen the expectation-value of the 
photon number to be the same as in the Fock case in Fig. ^. We see that for coherent 
states F{t) is mostly negative (solid curve in Fig. ^) while for Fock states F{t) is mostly 
positive (solid curve in Fig. ^. This illustrates the surprising fact that the correlation 
between the modes can be even more quantum mechanical for coherent states as compared 
to initial Fock states with the same average number of photons. 



6 Squeezing and Detuning 



In this section we will consider the influence of detuning on squeezing. The squeezing 
will be defined with respect to the phase-dependent quadrature amplitude Xg(t) (see e.g. 
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Figure 8: The correlation function F{t) = f (t) / y^ {uait)) {nb{t)) when k'^ = 1.5. 
The sohd cure shows F{t) when the system is initiahy a coherent state \a)a ^ 1/5);, 
with a = \/50, P = VW and the dotted curve when a = \/50, /3 = 0. 



Ref.iT 



Xg{t) = -^ {a{t)e''' + h{t)e-''' + h.c] 
v2 



(52) 



A number of schemes to make quadrature phase measurements have been discussed in 
the literature (see e.g. Ref.|]l4|). These schemes involve homodyning the signal field with 
a reference signal known as the local oscillator. In this case 9 is the phase of the local 
oscillator. The variables Xo(t) and X^/2(^) are canonical, i.e. they obey the commutation 
relation 

[^0, X^/2] = -2i . (53) 

The corresponding uncertainty product therefore is 



AXoAX,/2 > 1 



(54) 



We now consider squeezing in terms of the observable X0{t) for initial Fock states or 
coherent states. 
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6.1 The Fock State 

If the system initially is in the Fock state \r, s) the variance in Xe(t) is expressed by 



Var[Xeit)] = \Teit)mr + s + 1) , 



(55) 



where 



|T,(t)|2 = |e-^SW +^0- A.it) e-Mt) - ^^p 
= x{t) [1 + y{t) - 2 ( cos{nt - 2e)G(t) + sin(fit - 29)H{t) 
Here we have defined the functions 



(56) 



G{t) 



rni ;- tan((yft-\/A;^ — 1) [1 — tanh 5\ 

tan^ {giVk"^ - 1) + tanh^ 6 



r. rj tanh((7tVl — A;^) [1 + tan^ 7] 

1 + tanh^(5ftVl - P) tan^ 7 



if k^>\ 



if fc2 < 1 



(57) 



and 



Hit) 



r^rri \coi\i5 - tanh (5 [I +tw?{gt^k^ - I)] . 
^ tan2(c/tv/p^T)+tanh2 5 ^ 



/^— p .^ _ tan 7 [1-tanh (^tVl - F)] . 
^ 1 + tanh^ (^tVl - A;2) tan^ 7 ^ 



if A;2 > 1 



if fc2 < 1 



(5J 



x(t) and y{t) are given by Eq. (^) and Eq. (pif), respectively. The expressions for tan 7 
and coth 5 are, furthermore, given by Eq. ([T6|) and Eq. ([T8|) . The special case when fc^ = 
the variance in Eq. (E^) reduces to 



l/ar[Xo(t)] = e~'^*(r + s + 1) 



(59) 



and 



l/ar[X,/2(t)] = e''\r + s + 1) , (60) 

in agreement with well known results (see e.g. [0). We immediately see that this is a 
minimum-uncertainty state when r = s = 0. Therefore, changing the phase of the local 
oscillator, 6*, by 7r/2 enables one to move from enhanced to squeezed quadrature phase 
fiuctuations. 

By inspection of the actual expressions for A+(t) and A^{t) in Eq. (|l^ we see that 
Eq. ( |56D is periodic when k"^ satisfies Eq. (|^). The corresponding revival time is therefore 
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Figure 9: The thin sohd curve shows AXo(t), the dotted curve shows AXj^/2{t). 
The thick sohd curve is the corresponding uncertainty product. Ah curves corre- 
spond to the vacuum case and k"^ = 9/5. 



again given by Eq. (^Sf). Fig. |^ shows the uncertainties AXo(t) = ^yVar[XQ{t)] (thin 
solid curve) and AX-„/2{t) = A/V^ar]XjJy^(t)J (dotted curve) when the system is initially 
in vacuum for k"^ = 9/5. These curves oscillate with the period gUev = vry^ ~ 7.02, where 
we can choose p = 2 and n = 3. However, the period of the uncertainty product is just half 
the period of the individual uncertainties, i.e. uncertainty product oscillate with period 
gt = 7i\^/2 ^ 3.51. This is due to the fact that Var[XQ(t + trev/2)] = Var[XT^/2{t)], i.e. 
the curve Var[XT,/2(t)] is just shifted trevf^ relative to Var[XQ(t)]. 

Due to revivals the uncertainty product is one at the revival times. The uncertainties 
AXo(t) and AX^/2(^) are then equal and there is no squeezing. We now observe that we 
can write AXQ(t) = c(t) + d{t) and AX^/2('^) = c(t) —d(t), i.e. the uncertainty product is 
AXg(t)AX^,2(^) = c^(^) ~ d'^{t)j where the functions c(t) and d{t) can be obtained from 
Eq. (1561). From this observation we infer that extremal values of AXo(t) and AX7r/2(i) 
are associated with extremal values of the uncertainty product. This fact is also clearly 
illustrated in Fig. ^ 



Let us now study the case k < 1. Fig. [TO shows a logarithmic plot of the uncertainties 
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Figure 10: The thin sohd curve shows AXo(t), the dotted curve shows AX^/2(t). 
The thick sohd curve is the corresponding uncertainty product. All curves corre- 
spond to the vacuum case and A;^ = 0.5. 



AXo(t) and AX^/2(^) and the corresponding uncertainty product when the system is 
initially in the vacuum state and k^ = 0.5. As we can see, the uncertainties increase 
essentially at an exponentially rate as a function of time. However, the curves also exhibit 
an oscillatory behavior with various local minima. These minima actually occur when the 



time-derivative of Eq. (SB) is equal to zero. When gt^/1 — k? <^ 1 the derivative of 
this equation is complicated and the various minima must be determined numerically. 
When gtyjl — k"^ ^ 1 the periodicity for large gt is determined by the the cosine- and 
sine- functions in Eq. (p6D, i.e. the local minima of AXo(t) and AX7r/2(i) have a period 
gt = ir/k. The period of the uncertainty product is half of this period, as discussed above. 



6.2 The Coherent State 



The coherent case, when the a-and 6-mode are initially in the coherent state |a,/?). 
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the variance of X0{t) is given by 

Var[Xe{t)] = \Te{t)\\ (61) 

where again |Te(t)p is given by Eq. (|56|) . The most striking feature of Eq. (|6l|) is the 
fact that it is independent of the initial values of a and (3. Moreover, the only difference 
between the Fock variance in Eq. (^5]) and the coherent variance in Eq. (0) is the factor 
r + s + 1. The coherent state variances are therefore equal to the vacuum state variances. 



7 Signal-to-Noise Ratio 

We define the signal-to-noise ratio Pait) for the a- mode by 

where Uait) = a\t)a(t), see Eq. (^91) . A similar definition for the 6-mode holds. To get a 
good signal we want such ratios to be as large as possible. Let us suppose the system is 
initially in the Fock state \r, s). The ratio Pa{t) is then given by 

Pait) = r + Mt){r + s + l) 1 3 

^/no{t) + nl{t) V2rs + r + s + l 

where no{t) is given by Eq. (|35|) . The small and large time limits of Pa{t) are 



Pait) 



r + s 



^/2rs + r + s + l 
r + {gt)'^{r + s + l) 



if t ^ oo , A;^ < 1 , 

(64) 
if t ^ , for all k^ . 



, gty/2rs + r + s + 1 
If r = we then have a well-defined limit for small times, i.e. Pa(0) = (see e.g. 



Fig. |TT] (thick curves)). On the other hand, if r 7^ then Pa(0) = 00 (see e.g. Fig. |TT 
(thin curves)). This is due to the fact that a Fock state has no variance. The behavior of 
Pa{t) at t > depends strongly on k"^. li k^ < 1 the signal-to-noise ratio approaches the 
value as given in Eq. (l6^). If, on the other hand, fc^ > 1 the ratio Pa{t) will oscillate. The 
period of Pa{t) is determined by Eq. (|35D , i.e. the revival time is 
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Figure 11: The thin curves show Pa{t) when r = 1, s = 1. The thin sohd curve 
corresponds to k'^ = 1.5 and the thin dotted one corresponds to fc^ = 0.5. The 
thick curves show Pa{t) when r = 0, s = 10. The thick sohd curve corresponds to 
k'^ = 1.5 and the thick dotted one corresponds to k'^ = 0.5. The revival time is 
gtrev = V^ir w 4.44. 



gtrevvk 



niT 



(65) 



where n is an positive integer. The value of Pa{t) at the revival time is or oo, depending 
on the value of r as discussed above (see e.g. Fig. |Tl|). Therefore, we can achieve an 
infinite signal-to-noise ratio at the revival times by choosing r 7^ 0. 

Another interesting feature of the signal-to-noise ratio Pa{t) is that it can be approx- 



imatively constant between revival times, which also is illustrated in Fig. [Tl|. This is to 
be compared with the oscillatory behavior of e.g. pii(t) (piiit) in Fig. |l| is oscillating fast 
for times where Pa{t) is approximately constant). 

Let us now study the extrema of the signal-to-noise ratio, i.e. when the time derivative 
of Eq. (pB|) is zero. Pa{t) has always an extremum when 



gtVk^ - 1 = mi/2 



(66) 
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Figure 12: The thin curves show pa{t) when r = 100, s = 1. The thin sohd curve 
corresponds to k"^ = 1.5 with the global minimal value /o^^*'" = 152/v452 ~ 7.15. 
The thin dotted one corresponds to A;^ = 0.5. The thick curves correspond to the 
case r = 1, s = 100. The thick solid curve show pa{t) when k'^ = 1.5 with the local 
maximal value p^^'*'" = 205/(2\/452) « 4.82. The thick dotted one corresponds to 
k"^ = 0.5. Both dotted curves approach the value pa = 102/\/302 ^ 5.87. 



where n is an odd integer. The value Pait) at these particular times is 

rfc^ + s + 1 



extr 
Pa 



A/P(2rs + r + s + 1) 



If, in addition, r, s and k'^{> 1) satisfy the inequality 



< 



< 



1 



(67) 



(68) 



s — r + 1 /c^ — 1 

the extremum in Eq. (p7| ) is a local maximum. This inequality is e.g. satisfied when 
r = 1, s = 100 and A;^ = 1.5 as shown in Fig. |T2| (thick solid curve). The local maximum 
in this case is p^^*^ = 205/(2-\/452) ^ 4.82. Moreover, when Eq. ( |68D is satisfied the ratio 
Pa{t) has two other extrema as well. These extrema are global minima and occur when 



gWk"^ - 1 = ± 



arcsm 



(P 



+ 72^71 



(69) 
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where n_|_ = 0, 1, 2, 3, ... and n_ = 1, 2, 3, .... The particular value of this minimuiii is 



independent of the detuning parameter k'^. The minimum value corresponding to the 
thick solid curve in Fig. is p™'' = 2^/101/302 ^ 1.16. 

If, on the other hand, Eq. (|68|) is violated then p^^**" in Eq. ( |67D is the only extremum. 



It is a global maximum when r = (see e.g. thick solid curve in Fig. 11) and a global 
minimum when r 7^ (see e.g. thin solid curve in Fig. |1^). In both cases the global 
extremum is given by Eq. (|67D. 

When the system is initially in a coherent state the signal-to-noise ratio is essentially 
the same as in the Fock case. The essential difference is that Pa{t) is finite at the revival 
times since the variance is not zero. 

An alternative definition of the signal-to-noise ratio can be obtained by considering 
the quadrature operator 

Xa{t) = ^UJ[a + a^Ui , (71) 

where Uj is the time-evolution operator in Eq. (plQl) . The signal-to-noise ratio is then 
defined by 

Var[Xa{t)\ 
From this definition we immediately see that T]a{t) = for a Fock state. Eq. (|72D is, 
however, non-trivial for a coherent state \a)a <S) \P)b = |a,/3). In this case the ratio ria{t) 

is 

^ K{t) + x{t) { \a\' - 2Re[a(3A.{t)] + \(3\Mt) } .^3^ 

"^"^ ' no(t) + 1/2 ' ^ ' 

where 

K{t) = Re [e-2^S(*)(a _ f3*A*_{t)y] . (74) 

The functions x{t), y(t) and no(t) are given by Eq. (|33|) , Eq. (|3^) and Eq. (|35|) , respec- 
tively. The signal-to-noise ratio in Eq. (|73|) cannot take on any value: it has an upper 
bound. According to Yuen PB[ the maximum signal-to-noise ratio obtainable for the 



(two-photon) coherent state |a,/5) cannot exceed the value 4(na(t)) [{na{t)) + 1], i.e. 

Va{t) < 4(n„(t)) [{na{t)) + 1] , (75) 

where {na{t)) = {(3 , a\a'^ {t)a{t)\a, f3) . It is natural to ask for what values of a, (3 and k"^ 
can we reach this optimal Yuen- limit ? As an illustrative example we choose a large /c^, 
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Figure 13: The signal-to-noise ratio ria{t) when k'^ = 10. The system is initially 
in the coherent state |a,/3) where a = and /3 = 3 (solid curve). The dotted curve 
is the corresponding Yuen-limit. 

e.g. k"^ = 10, and a = and (3 = 3. The result of a numerical calculation is shown in 
Fig. |TB|. As we can see from such a calculation, there is approximately a factor of 2.5 
in difference between the maximum ria{t)-value (solid curve) and the optimal Yuen-limit 
(dotted curve). In other words, the signal is not very far from its optimal value. The 
price we have to pay to come this close to the optimal Yuen-limit is the low value of ria{t) 
itself. If we, on the other hand, choose a small k"^, e.g. k"^ = 1.5 (with the same a and (3), 
the physics is different. The maximum ratio T]a{t) is then (9(10^), which means that the 
signal is strong as compared to the noise. In this case the Yuen-limit is, however, (9(10^), 
so rjait) is very far from the optimal Yuen- value. Numerical studies suggest that these 
observations are quite generic. 

8 Final Remarks 

By means of an, in general time- dependent, Bogoliubov transformation or what is 
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called a two-mode squeeze operator S{z) in quantum optics |]TT|, given by 



S{z) = exp [z*ab — za^b^) , (76) 

where z = rexp{i(f)), one can diagonalize the Hamiltonian Eq. (Qj instantaneously in 
terms of the canonical operators 

A(z) = S{z)aS\z) = acoshr + 6' exp(20) sinhr , 

B{z) = S{z)bS\z) = fecoshr + a^exp(i0) sinhr . (77) 

With the choice Eq. (|T2|) we e.g. find that 

H = nAA^{z)A{z) + nBB^{z)B{z) + fio (78) 

where 



Qa = uj-+uj+\ I- g^/ujl , 



Hb = —UJ- + UJ- 



Vi-w^' 



fio = u+^l-g^/ul , (79) 

and where uj± = {ua ± ujb)/^. The squeeze parameters are given by 

cosh(2r) = l/y'T^TM , = | - ^^ • (§0) 

The diagonalization procedure above is well-defined as long as g"^ < 00"^. If, on the other 
hand, g^ > u"^ then the instantaneous eigenvalues Qa and Qb become complex. This 
indicates an instability in the system. As we have seen in the present paper, one can 
explicitly solve the time-dependent Schrodinger equation in such a situation and, as we 
have seen above, the physics of the instability actually manifest itself in an exponentially 
increasing particle, i.e. photon, production. In the case /c^ < 1 we have also seen that the 
asymptotic time-limit of the reduced systems can be described by thermal distributions. In 
general the distributions are not identical (see e.g. Eqs.(p9D and (^0])). But, nevertheless. 



if the initial state is a pure state, the reduced entropies are the same (see Refs. [^ ^ |35| ) 
and the entropy of the whole system is, of course, zero. 

Some features of the simple model of parametric down-conversion as studied in the 
present paper also enters in models considered in quantum cosmology as e.g. the phenom- 



ena of non-adiabatic transitions (see e.g. Ref.[^]). To the extent that such an analogy 
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holds we can therefore simulate some aspects of such models in the laboratory by means 
of non-linear quantum optics. 
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